Mixed and discontinuous finite volume element schemes for the optimal control of immiscible flow in porous media

We introduce a family of hybrid discretisations for the numerical approximation of optimal control problems governed by the equations of immiscible displacement in porous media. The proposed schemes are based on mixed and discontinuous finite volume element methods in combination with the optimise-then-discretise approach for the approximation of the optimal control problem, leading to nonsymmetric algebraic systems, and employing minimum regularity requirements. Estimates for the error (between a local reference solution of the infinite dimensional optimal control problem and its hybrid approximation) measured in suitable norms are derived, showing optimal orders of convergence

Competing mechanisms of stress-assisted diffusivity and stretch-activated currents in cardiac electromechanics

We numerically investigate the role of mechanical stress in modifying the conductivity properties of the cardiac tissue and its impact in computational models for cardiac electromechanics. We follow a theoretical framework recently proposed in [Cherubini, Filippi, Gizzi, Ruiz-Baier, JTB 2017], in the context of general reaction-diffusion-mechanics systems using multiphysics continuum mechanics and finite elasticity. In the present study, the adapted models are compared against preliminary experimental data of pig right ventricle fluorescence optical mapping. These data contribute to the characterization of the observed inhomogeneity and anisotropy properties that result from mechanical deformation. Our novel approach simultaneously incorporates two mechanisms for mechano-electric feedback (MEF): stretch-activated currents (SAC) and stress-assisted diffusion (SAD); and we also identify their influence into the nonlinear spatiotemporal dynamics. It is found that i) only specific combinations of the two MEF effects allow proper conduction velocity measurement; ii) expected heterogeneities and anisotropies are obtained via the novel stress-assisted diffusion mechanisms; iii) spiral wave meandering and drifting is highly mediated by the applied mechanical loading. We provide an analysis of the intrinsic structure of the nonlinear coupling using computational tests, conducted using a finite element method. In particular, we compare static and dynamic deformation regimes in the onset of cardiac arrhythmias and address other potential biomedical applications

Stabilized mixed approximation of axisymmetric Brinkman flows

This paper is devoted to the numerical analysis of an augmented finite element approximation of the axisymmetric Brinkman equations. Stabilization of the variational formulation is achieved by adding suitable Galerkin least-squares terms, allowing us to transform the original problem into a formulation better suited for performing its stability analysis. The sought quantities (here velocity, vorticity, and pressure) are approximated by Raviart−Thomas elements of arbitrary order k ≥ 0, piecewise continuous polynomials of degree k + 1, and piecewise polynomials of degree k, respectively. The well-posedness of the resulting continuous and discrete variational problems is rigorously derived by virtue of the classical Babuška–Brezzi theory. We further establish a priori error estimates in the natural norms, and we provide a few numerical tests illustrating the behavior of the proposed augmented scheme and confirming our theoretical findings regarding optimal convergence of the approximate solutions

Finite element and finite volume-element simulation of pseudo-ECGs and cardiac alternans

In this paper, we are interested in the spatio-temporal dynamics of the transmembrane potential in paced isotropic and anisotropic cardiac tissues. In particular, we observe a specific precursor of cardiac arrhythmias that is the presence of alternans in the action potential duration. The underlying mathematical model consists of a reaction–diffusion system describing the propagation of the electric potential and the nonlinear interaction with ionic gating variables. Either conforming piecewise continuous finite elements or a finite volume-element scheme are employed for the spatial discretization of all fields, whereas operator splitting strategies of first and second order are used for the time integration. We also describe an efficient mechanism to compute pseudo-ECG signals, and we analyze restitution curves and alternans patterns for physiological and pathological cardiac rhythms

Nonlinear diffusion & thermo-electric coupling in a two-variable model of cardiac action potential

This work reports the results of the theoretical investigation of nonlinear dynamics and spiral wave breakup in a generalized two-variable model of cardiac action potential accounting for thermo-electric coupling and diffusion nonlinearities. As customary in excitable media, the common Q10 and Moore factors are used to describe thermo-electric feedback in a 10-degrees range. Motivated by the porous nature of the cardiac tissue, in this study we also propose a nonlinear Fickian flux formulated by Taylor expanding the voltage dependent diffusion coefficient up to quadratic terms. A fine tuning of the diffusive parameters is performed a priori to match the conduction velocity of the equivalent cable model. The resulting combined effects are then studied by numerically simulating different stimulation protocols on a one-dimensional cable. Model features are compared in terms of action potential morphology, restitution curves, frequency spectra and spatio-temporal phase differences. Two-dimensional long-run simulations are finally performed to characterize spiral breakup during sustained fibrillation at different thermal states. Temperature and nonlinear diffusion effects are found to impact the repolarization phase of the action potential wave with non-monotone patterns and to increase the propensity of arrhythmogenesis

Temperature dependence of the anomalous effective action of fermions in two and four dimensions

The temperature dependence of the anomalous sector of the effective action of fermions coupled to external gauge and pseudo-scalar fields is computed at leading order in an expansion in the number of Lorentz indices in two and four dimensions. The calculation preserves chiral symmetry and confirms that a temperature dependence is compatible with axial anomaly saturation. The result checks soft-pions theorems at zero temperature as well as recent results in the literature for the pionic decay amplitude into static photons in the chirally symmetric phase. The case of chiral fermions is also considered.Comment: RevTex, 19 pages, no figures. References adde

Anomalous mesonic interactions near a chiral phase transition

Using constituent quarks coupled to a linear sigma model at nonzero temperature, I show that many anomalous mesonic amplitudes, such as $\pi^0 \rightarrow \gamma \gamma$, vanish in a chirally symmetric phase. Processes which are allowed, such as $\pi^0 \sigma \rightarrow \gamma \gamma$, are computed to leading order in a loop expansion.Comment: 4 pages, REVTeX, submitted to Phys. Rev. Let

How pi0 -> gamma gamma changes with temperature

At zero temperature, in the chiral limit the amplitude for pi0 to decay into two photons is directly related to the coefficient of the axial anomaly. At any nonzero temperature, this direct relationship is lost: while the coefficient of the axial anomaly is independent of temperature, in a thermal bath the anomalous Ward identities do not uniquely constrain the amplitude for pi0 -> gamma gamma. Explicit calculation shows that to lowest order about zero temperature, this amplitude decreases.Comment: 12 pages, ReVTeX, 4 figures, to be published in Phys. Rev. D. New section 5 with proof of the Adler-Bardeen theorem at low

Measurement of the relative rate of prompt χc0, χc1 and χc2 production at √s=7TeV

Prompt production of charmonium χc0, χc1 and χc2 mesons is studied using proton-proton collisions at the LHC at a centre-of-mass energy of √s=7TeV. The χc mesons are identified through their decay to J/ψγ, with J/ψ→μ+mu− using photons that converted in the detector. A data sample, corresponding to an integrated luminosity of 1.0fb−1 collected by the LHCb detector, is used to measure the relative prompt production rate of χc1 and χc2 in the rapidity range 2.0<y<4.5 as a function of the J/ψ transverse momentum from 3 to 20 GeV/c. First evidence for χc0 meson production at a hadron collider is also presented

Measurement of χ c1 and χ c2 production with s√ = 7 TeV pp collisions at ATLAS

The prompt and non-prompt production cross-sections for the χ c1 and χ c2 charmonium states are measured in pp collisions at s√ = 7 TeV with the ATLAS detector at the LHC using 4.5 fb−1 of integrated luminosity. The χ c states are reconstructed through the radiative decay χ c → J/ψγ (with J/ψ → μ + μ −) where photons are reconstructed from γ → e + e − conversions. The production rate of the χ c2 state relative to the χ c1 state is measured for prompt and non-prompt χ c as a function of J/ψ transverse momentum. The prompt χ c cross-sections are combined with existing measurements of prompt J/ψ production to derive the fraction of prompt J/ψ produced in feed-down from χ c decays. The fractions of χ c1 and χ c2 produced in b-hadron decays are also measured